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ABSTRACT 

We study the metallicities and abundance ratios of early-type galaxies in cosmologi- 
cal semi-analytic models (SAMs) within the hierarchical galaxy formation paradigm. 
To achieve this we implemented a detailed galactic chemical evolution (GCE) model 
and can now predict abundances of individual elements for the galaxies in the semi- 
analytic simulations. This is the first time a SAM with feedback from Active Galactic 
Nuclei (AGN) has included a chemical evolution prescription that relaxes the instan- 
taneous recycling approximation. We find that the new models are able to reproduce 
the observed mass-metallicity (Mi,-[Z/H]) relation and, for the first time in a SAM, 
we reproduce the observed positive slope of the mass- abundance ratio (Af^-[a/Fe]) 
relation. Our results indicate that in order to simultaneously match these observa- 
tions of early-type galaxies, the use of both a very mildly top-heavy IMF (i.e., with a 
slope of a: = 1.15 as opposed to a standard x = 1.3), and a lower fraction of binaries 
that explode as Type la supernovae appears to be required. We also examine the rate 
of supernova explosions in the simulated galaxies. In early-type (non-star forming) 
galaxies, our predictions are also consistent with the observed SNe rates. However, in 
star-forming galaxies, a higher fraction of SN la binaries than in our preferred model 
is required to match the data. If, however, we deviate from the classical model and 
introduce a population of SNe la with very short delay times, our models simultane- 
ously produce a good match to the observed metallicities, abundance ratios and SN 
rates. 

Key words: galaxies: formation - galaxies: chemical evolution 



1 INTRODUCTION 

The chemical properties and abundance ratios of galax- 
ies provide important constraints on their formation 
histories. Galactic chemical evolution has been mod- 
elled in detail in the monoli th ic collapse s c enario 
(e.g.. iMatteucci fc Greggiol Il986l: iFrancois et al.l |2004 



iRomano et al.ll2005l : |Pipino fc Matteuccill2004l2006l '). These 
models have successfully described the abundance distri- 
butions in our Galaxy and other spiral discs, as well as 
the trends of metallicity and abundance ratios observed in 
early-type galaxies. In the last three decades, however, the 
paradigm of hierarchical assembly in a Cold Dark Matter 
(CDM) cosmology has revised the picture of how structure 
in the Universe forms and evolves. In this scenario, galax- 
ies form when gas radiatively cools and condenses inside 
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dark matter haloes, w hich themselves follow dissipationless 
gravi tational collapse (|White fc Reeslll978l : IWhite fc FrenkI 
[199J). The CDM picture has been successful at predicting 
many observed properties of galaxies, though many poten- 
tial problems and open questions remain. It is therefore in- 
teresting to see whether chemical evolution models, when 
implemented within this modern cosmological context, are 
able to correctly predict the observed chemical properties of 
galaxies. 

The semi-analytic approach provides a cosmological 
framework in which to study galaxy formation and chemical 
evolution in different environments, by following the merger 
history of dark matter haloes and the relevant physical pro- 
cesses such as gas coo ling, star for mation and fe e dback (e.g. 



White fc Frenk 



I991I: fK auffma nn et all Il993l: ICole et all 



1994 12OOOI: iHatton et aTll2003i: iSoinerville fc Primacklll999l: 



Somcrvill e~et all bOOll ). A major challenge for models of 
galaxy formation within the CDM picture arises from the 
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mismatch between the shape of the mass function of the 
dark matter haloes and that of the baryonic condensations 
that we c all galaxies ("White & Frcnk 1991; Kauffma nn et alJ 
ll993l : [Somerville fc Primack 1999 : Benson c t al. 200^. The 
CDM theory predicts a steeper slope for low-mass halos, 
and a more gradual drop-off in the abundance of high-mass 
halos than is seen in luminous galaxies, implying that the 
formation of stars must be ineffi cient in both low-mass and 
high- mass haloes (|Moster et al.|[2009i '). However, the inclu- 
sion of physically motivated, if still ad hoc, feedback pro- 
cesses in the semi-analytic models can cure these discrep- 
ancies. The faint end of the luminosity function can be 
matched with a combination of supernova feedback and sup- 
pression of gas cooling in low mass haloes as a result of a 
photo-ionising background. At the bright end, heating by 
giant radio jets powered by accreting black holes has be- 
come a favored mechanism for preventing ove r-cooling and 
quen c hing star formation in massive halos (ICroton et al.l 
l2006l : iBower et al] l2006l : ISomerville et al.l |2008| ). This lat- 
est generation of semi-analytic models ('SAMs') is suc- 
cessful at reproducing many properties of galaxies at the 
present and at high redshift, for example, the luminosity 
and stellar mass function of galaxies, color-magnitude or 
star formation rate vs. stellar mass distributions, relative 
numbers of early and late-type galaxies, gas fractions and 
size distributions of spiral galaxi es, and the global star for- 
mation history (e.g . ICrot on et al. 2006: Bower ct al. 2006: 



Cattaneo et al.l 2006: Do L ucia et al.li 2006: So merville et al.l 
20081 : iKimm et al.ll2009i : iFontanot et al.ii2009i . to name just 



a few). 

The modelling of chemical enrichment of the galaxies 
and intergalactic (and intracluster) gas, however, has not 
been thoroughly developed in semi-analytic models, and to 
date most SAMs have only used the instantaneous recycling 
approximation (in essence only considering enrichment by 
type II supernovae) and trace only the total metal con- 
tent. There are, however, a few models that have in cluded a 
more refined treatment of the chemical enrichment. iThomad 
(|l999l) and iThomas fc KauffmannI (|l999l ') were the first to 
include enrichment by type la supernovae (SNe la) in mod- 
els with cosmologically motivated star formation histories. 
However, rather than implementing the chemical evolution 
self-consistently within a semi-analytic model, they made 
use of star formation histories from the SAM and assumed 
a closed-box model (no gas inflows or outflows) for the chem- 
ical evolution. They calculated the evolution of [Fe/H] and 
[Mg/Fe] and found a decreasing trend of [Mg/Fe] with in- 
creasing galaxy luminosity, in stark disagreement with o b- 



servations 
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, .Worthev et all Il992l : iTrager et af] l2000al lbl: 
20051 ). 



The first semi-analytic model to self-consistently track a 
variety of elements due to enric hment by SNe la and type II 
supernovae (SNe II) was that of lNagashima et a l. (2005a b). 
Among other things, they adopted a bimodal IMF described 
by a standard IMF for normal quiescent star formation in 
discs and an extremely flat 'top-heavy' IMF during merger- 
driven starbursts. This model was motivated by the diffi- 
culty that semi-analytic models with a standard IMF expe- 
rienced in reproducing the observed pop ulation of very lum i- 
nous sub-mm galaxies at high redshift (|Baugh et al Tl2005h . 
However, the notion that early-type galaxies form their stars 
with an IMF flatter than standard is not new and has been 



proposed many times in the past as a plausible explanation 
for the abundance patterns in early-type galaxies and in the 
ICM of galaxy clusters (e.g., Worthey, Faber & Gonzale; 
19921: iMatteucci fc GibsoDlll995..: Gibson fc Mattcucci 199 



Thomas. Greggio fc Bended 11999^ 7 The predictions of the 
Nagashima et al. model were in good agreement with the 
abundances of the intracluster medium (ICM) of galaxy clus- 
ters, matching the trend of individual elements (O, Fe, Mg, 
Si) and abundance ratios with ICM temperature. However, 
the same model failed to reproduce the trend of [a/Fe] in 
early-type galaxies, where they found that the abundance 
ratio decreases with increasing galactic velocity dispersion, 
again in clear contradiction with observations. Very recently, 
IPipino et all jioOS) have coupled gala ctic chemical evolu - 
tion to the GallCS semi-analytic model llHatton et al.ll2003l), 
and obtained results similar to those of Nagashima et al] 
(2005b). 

In a simple closed-box picture, it is well-known that 
galaxies with short star formation timescales are expected 
to have enhanced [a/Fe] ratios (because their enrich- 
ment is dominated by a-rich Type II SNe), while galax- 
ies with extende d star formation histories tend to have 
lower [g/ Fe] (e . g. IWorthev et alj|l992l; iThornas et alj|l999l: 
iThomasI 1 19991 : iThomas fc KauffmannI 1 19991 : iTrager et all 
2000bj) , because of the additional Fe contributed by de- 
layed Type la enrichment. Therefore a possible interpreta- 
tion of the difficulties that CDM-based galaxy formation 
models have experienced in reproducing the positive trend 
between mass, luminosity, or velocity dispersion and abun- 
dance ratio is related to the issue of so-called "downsiz- 
ing". This refers to the variety of observational evidence 
that high-mass galaxies formed their stars early and over 
short timescales, while low-mass ga laxies have more ex- 
tended star formation histories (see iFontanot et alj |2009| . 
for a summary). Before the inclusion of AGN feedback 
or some other mechanism that quenches star formation in 
massive halos, CDM-based galaxy formation models pre- 
dicted the opposite trend (massive galaxies continued to 
accrete gas and form stars until the present day, leading 
to extended star formation histories). It has been demon- 
strated that including radio-mode AGN feedback in semi- 
analytic models leads to a "downsizing" trend for star for- 
mation that is at least qualitatively i n better agreement with 
observations ( De Lucia et al. 2006; Somerville et al. 2008|; 
ITrager fc Somerville. 2009 : .Fontanot ct al. 2009). Therefore 
we expect that the new models might do better at repro- 
ducing the trend of [a/Fe] with mass as the more massive 
galaxies will have shorter star formation timescales. 

Clearly, observations of chemical abundances and abun- 
dance ratios in various phases (stellar, ISM, ICM) offer the 
opportunity to obtain strong constraints on galaxy forma- 
tion histories and the physics that shapes them. However, 
in order to take advantage of these observations, it is neces- 
sary to implement detailed modeling of chemical evolution 
into a full modern SAM that includes the relevant phys- 
ical processes (e.g. triggered star formation and morpho- 
logical transformation of galaxies via mergers, the growth 
of supermassive black holes, and AGN feedback). In this 
work we incorporate detailed chemical e volution into the 
semi-analytic galaxy formation model of ISomerville et al] 
(2008), taking into account enrichment by SNe la, SNe II 
and long-lived stars, and abandoning the instantaneous re- 
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cycling approximation by considering the finite lifetimes of 
stars of all masses. The delay in the metal enrichment by 
SNe la is calculated self-consistently according to the life- 
times of the progenitor stars. This is, to our knowledge, the 
first time that detailed chemical evolution has been included 
in a semi-analytic model with AGN feedback (both radio- 
mode heating and AGN-driven winds). The base model in- 
cludes gas inflows due to radiative cooling of gas and out- 
flows due to supernova and AGN-driven winds. We compute 
the abundances of many a and Fe-peak elements for early- 
type galaxies of different masses, exploring different IMF 
slopes and values for the fraction of binaries that yield a 
SN la event, and compare these with observations of abun- 
dances and abundance ratios for a sample of local early-type 
gal axies. We also calcu l ate S Ne rates using both the classi- 
cal iGreggio fc Renzinil (Il983 l) approach for type la SN and 
the more recen t Delay-Time-Distribution (DTD) formalism 
l|Greggidl2005l ) . Another improvement in the present work is 
our use of re-calibrated estimates for chemical abundances 
obtained from line-strengths in early type galaxies (see Ap- 
pendix B for details). 

The outline of the paper is as follows. In Section 2 we 
give an overview of the main ingredients of the semi-analytic 
model. In Section 3 we describe in detail the adopted treat- 
ment for the chemical evolution. In Section 4 we present our 
predictions and compare them with observations. In Section 
5 we summarise our flndings and present our conclusions. 
Two appendices describe the detailed implementation of the 
chemical evolution model and the data used in this paper. 



2 THE SEMI-ANALYTIC MODEL 

In this section we summarise the basic ingredients of the 
SAM used to model the formation and evolution of galax- 
ies. These include the growth of structure of the dark 
matter component in a hierarchical clustering framework, 
radiative cooling of gas, star formation, supernova feed- 
back, AGN feedback, galaxy merging within dark matter 
haloes, metal enrichment of the ISM and ICM, and the 
evolution of stellar populati o ns. The reader is referr ed to 
ISomerville fc Primackl ll 19991'). ISomerviUe et all (|200ll ') and 
especiallv ISomerville et al.l (|2008l , hereafter SOS) for a com- 
prehensive and detailed description of the different prescrip- 
tions used in this semi-analytic model. In what follows we 
briefly sketch the modelling of the most important physical 
processes. 



2.1 Dark matter merger trees and galaxy merging 

The merging histories (or merger trees) of dark mat- 
ter haloes are constructed based on the Extended Press- 
Schechter formalism usin g the method described in 
ISomerville fc KolattI l|l999l ). with improvements described 
in SOS. Each branch in the tree represents a merger event, 
and, in order to make the process finite, the trees are fol- 
lowed down to a minimum progenitor mass of 10^° -^Q- 

Whenever dark matter haloes merge, the central galaxy 
of the largest progenitor becomes the new central galaxy, 
and all others become 'satellites'. Satellite galaxies may 
eventually merge with the central galaxy due to dynam- 
ical friction. To model the timescale of the merger pro- 



cess we use a variant of th e Chandrasekhar formula from 
iBovl an-Kolchin et al] l|2008l ). Tidal stripping and destruc- 
tion of satellites are also included as described in SOS. 



2.2 Gas cooling, star formation and supernova 
feedback 

Before the Universe is reionised, each halo contains a mass of 
hot gas equal to the universal baryon fraction times the virial 
mass of the halo. After reionisation, the photo-ionising back- 
ground can suppress the collaps e of ga s int o low-mass ha- 
los. W e use the results of lGnedinl (|200G| ) and lKravtsov et al.l 
(|2004 ) to model the fraction of baryons that can collapse 
into haloes of a given mass after reionisation. 

When a dark matter halo collapses, or merges with a 
larger halo, the gas within it is shock-heated to the virial 
temperature of the halo, and gradually radiates and cools 
at a rate given by the cooling function. To calculate this 
function we use the metallicity-dcpondent radiative cooling 
curves of .Sutherland fc Dopita, (199S ) . A detailed descrip- 
tion of how the cooling process is modelled can be found in 
SOS. The rate at which gas can cool is given by: 

1 ^cool 1 / \ 

nicooi = -Jnhot , (1) 

Z 7'vir f'cool 

where rrihot is the mass of the hot halo gas, Tvir is the virial 
radius of the dark matter halo, rcooi is the radius within 
which all of the gas can cool in a time tcooi , which itself de- 
pends on density, metallicity and temperature. In our mod- 
els, we assume that the cold gas is accreted only by the 
central galaxy of the halo, but in reality satellite galaxies 
should also receive some measure of new cold gas. This as- 
pect of the modelling should be improved (cf. iPipino et al.l 
[200^ and for the present study we restrict our analysis to 
only the central galaxy of each halo, except when otherwise 
stated. 

When the gas cools we assume that it settles into 
a rotationally supported disc. The radial sizes of the 
discs are calcula t ed ac cording to the results described in 
ISomerville et all (120081 ). and agree well with observed disc 
sizes to 2 ^ 2. 

We model the star formation rate in quiescent discs 
with a recipe bas ed on the empirical Schmidt-Kennicutt law 
(|Kennicuttlll9Sg| '): 

27rAxS^'«r^,, 



exp ( — Nk T-crit / r-gas ) 



(2) 



where rii* is the star formation rate, Eq 



a 



the average surface density of the cold gas, rgas is the scale- 
length of the gaseous disc (assumed to be an exponential 
disc with its scale-length proportional to that of the stellar 
disc) , rcrit is the radius at which the gas reaches the critical 
surface density threshold for star formation (Scrit), and Ak 
and Nk are the normalisation and slope of the SFR law. 
We adopt the values Ak = 8.35 x 10"'^, Nk = 1.4 and 
E„it = 6M0pc-^ as in SOS. 

Galaxy mergers in the SAM trigger enhanced episodes 
of star formation. The burst is modelled by two parameters, 
the time-scale and the efficiency of the burst. The time-scale 
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is a function of the virial velocity of the progenitor galax- 
ies, the equation of state of t he gas, the cold gas fr action 
in the discs, and the redshift ijRobertson et al]l2006t l. The 
efficiency, which is defined as the fraction of the cold gas 
reservoir (of both galaxies) that is turned into stars dur- 
ing the burst, is assumed to be a power-law function of the 
mass ratio of the merging galaxies, and the expo nent of the 
powe r-law depends on the galaxy morphology (|Cox et al.l 
I2OO8I ). The coUisional starburst occurs in addition to any 
ongoing 'normal' quiescent star formation, which continues 
uninterrupted through the merger but is usually insignifi- 
cant in comparison to the burst. Any new stars formed dur- 
ing the burst are always placed in the bulge component of 
the resulting galaxy. 

As supernovae occur, they inject energy into the ISM 
and reheat the cold gas, which is then expelled from the 
disc and incorporated into the hot halo gas where it can 
cool again. The rate of reheating by SNe is given by 

SN ( Vdisc \ . 

where and arh are free parameters. The circular velocity 
of the disc Vdisc is taken to be equal to the maximum rota- 
tional velocity of the dark matter halo. Some fraction of the 
reheated gas can also be ejected from the halo entirely into 
the diffuse Intergalactic Medium. This fraction is described 
by: 

/cjcct(Kir) = [1.0 -I- (Vvir/Kjcct)""^""^']"', (4) 

where Qcjcct = 6 and Vcjoct is a free parameter in the range 
~ 100 — 150 km/s. This ejected gas is allowed to re-collapse 
into the halo at later times and once again becomes available 
for cooling. 



halo is seeded with a ~ 100 M© black hole, and these black 
holes are able to grow via two different accretion modes. 
The first accretion mode is fuelled by cold gas that is driven 
into the nucleus of the galaxy by mergers. This mode is ra- 
diatively efficient, and the accretion rates are close to the 
Eddington limit. Because this accretion mode is associated 
with optically bright classical quasars and AGN, it is re- 
ferred to as 'bright mode' or 'quasar mode' accretion. The 
second mode is fuelled by hot gas in a quasi-hydrostatic halo, 
and the accretion rate is modelled via the Bondi-Hoyle ap- 
proximation. Accretion rates in this mode are significantly 
sub-Eddington (~ 10~* to 10""^ times the Eddington rate) 
and the accretion is assumed to be radiatively inefficient. 
This mode is, however, associated with the production of 
giant radio jets, and is therefore referred to as the 'radio 
mode'. 

Energy released during 'bright mode' activity can cou- 
ple with the cold gas in the galaxy via radiation pressure, 
driving galactic scale winds that can eject cold gas from the 
galaxy. The mass outfiow rate due to the AGN driven wind 
is modelled by the following formula: 

^^agn — f^wind^^rad TT ^^acc , (5) 

''CSC 

where twind is the effective coupling efficiency, Ksc is the 
escape velocity of the galaxy and rhacc is the accretion rate 
of mass onto the black hole. 

The radio jets produced by 'radio mode' activity are as- 
sumed to inject thermal energy into the hot halo gas, partly 
or completely offsetting the cooling fiow. This process is re- 
sponsible for quenching the star formation in massive galax- 
ies (which contain massive black holes) and solves the 'over- 
cooling problem' that plagued CDM-based galaxy formation 
models for many years. 



2.3 Formation of Spheroids 

In most semi-analytic models each merger is classified as 
'major' or 'minor' depending on whether the ratio of the 
smaller to the larger galaxies' baryonic masses is greater 
than or less than the parameter /dup ~ 0.25, respectively. 
The usual assumption is then that, in a major merger, 
the bulge and disc stars of both progenitor galaxies, as 
well as the stars formed in the merger driven starburst 
(see below), are transferred to the bulge component of the 
resulting galaxy. In a minor merger, all the pre-existing 
stars of the smaller galaxy end up in the disc of the post- 
merger galaxy, and all the newly formed stars are placed 
in bulge. We follow a similar practise here, but instead of 
using a sharp threshold to define major or minor merg- 
ers, we use a more gradual transition function. In detail, 
when two galaxies with bulge masses B\ and -62, and disc 
masses D\ and D2 merge, the resulting galaxy has a bulge 
mass -Bncw ~ B\ -\- B2 -\- /sph(^i + B)2) and a disc mass 
I^ncw = (1 — /sph)(£'i -I- D-2). The value /sph is a continuous 
function of the total mass ratio (baryons and dark matter) 
in the central parts of the galaxy (see SOS). 

2.4 Black Hole Growth and AGN Feedback 

The models of SOS also track the growth of super-massive 
black holes and the energy they release. Each top-level DM 



2.5 Stellar Population Synthesis and Dust 

In order to compare the luminosities and colours of the 
galaxies in the simulations with real observations, we con- 
volve the star formation and chemical enrichment history of 
each galaxy with the m ulti-metallicity simple stel lar popu- 
lation (SSP) models of iBruzual fc Charloti (l2003h. We use 
the models based on the P adoval994 l|Bertelh et al.lll994 ) 
isochrones with a lChabriei] (|2001i ') IMF. 

We also model the effec t s of d ust extinction. Based on 
the model of ICharlot fc Falll (j200d ). we consider extinction 
due to two components, one due to the diffuse dust in the 
disc and another associated with the dense 'birth clouds' 
surrounding young star forming regions. The V-band, face- 
on extinction optical depth of the diffuse dust is given by 

TKO tX rdust,0-^cold?Tlcold/(''gas)^, (6) 

where rdust.o is a free parameter, Z^o\a is the metallicity of 
the cold gas, mcoid is the mass of the cold gas in the disc, 
and rgas is the radius of the cold gas disc. To compute the 
actual extinction we assign each galaxy a random inclination 
and use a standard 'slab' model. Additionally, stars younger 
than lO'^ yr are enshrouded in a cloud of dust with optical 
depth TBC,v = MBCTV.o, where /xbc = 3. Finally, to extend 
the extinction correction to o ther wavebands, we assume 
a Galactic attenuation curve (ICardelU et al.lll9S9l ) for the 
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diffuse dust component and a power-law extinction curve 
Ax (X (A/55OOA)", with n = 0.7, for the birth clouds. 



histories provided by the SAM. In the next subsection we 
describe the implementation of the new chemical evolution 
model in detail. 



2.6 Cosmological and Galaxy Formation 
Parameters 

We adopt a flat ACDM cosmology with = 0.2383, = 
0.7617, h = iyo/(100km s^^Mpc-i) = 0.732, ag = 0.761, 
and a cosmic baryon fraction of fb = 0.1746, following the re- 
sults of lSpergc l ct al. (2007). We adopt these parameters for 
consistency with the published models of SOS, but find that 
we obtain nearly identical results w ith the updated value s 
of the cosmological parameters from lKomatsu et"aLl l|2009t ). 

We leave the values of the free parameters associated 
with the galaxy formation models fixed to the fiducial values 
given in SOS. These values were chosen by requiring that the 
models reproduced key observations of nearby galaxies, such 
as the z ~ stellar mass function, and gas fractions and star 
formation rates as a function of stellar mass. These models 
have also been shown to produce reasonable ag reement with 
observ ed local galaxy colour distributions in iKimm et al.l 
l|2009f l. and with observed stellar mass functions and star 
form ation rates at high redshift (0 < z < 4; iFontanot et al.l 
I2OO9I ). In ij 14.11 we check that our new models, with the up- 
dated treatment of chemical evolution modelling, still repro- 
duce the key observational quantities with the same values 
of the free parameters. 



3 GALACTIC CHEMICAL EVOLUTION 

In SOS, the production of metals was tracked using a sim- 
ple approac h commonly ad opted i n semi-analytic mod- 
els (see e.g. ISomeryille fc P rimack 19991 : ICole eit all I2OOOI: 
iDe Lucia et al.l2004l : iKang et al.ii2005l ). In a given time-step, 
when we create a parcel of new stars dm*, we also create a 
mass of metals AMz ~ ydm,, which we assume to be in- 
stantaneously mixed with the cold gas in the disc. The yield 
y is assumed to be constant, and is treated as a free param- 
eterQ We track the mean metallicity of the cold gas Zcom, 
and when we create a new parcel of stars they are assumed 
to have the same metallicity as the mean metallicity of the 
cold gas in that time-step. Supernova feedback ejects metals 
from the disc, along with cold gas. These metals are either 
mixed with the hot gas in the halo, or ejected from the halo 
into the 'diffuse' Intergalactic Medium (IGM), in the same 
proportion as the reheated cold gas. The ejected metals in 
the 'diffuse gas' reservoir are also re-accreted into the halo 
in the same manner as the gas. 

In the present study, we discard the instantaneous re- 
cycling approximation and allow the ISM to be enriched by 
the products of type la and type II supernovae on their own 
timescales. Consequently, we now track individual elements, 
and not just the total metal content. The integrated ejecta 
of each element is not a free parameter, but instead is calcu- 
lated according to theoretical yields and the star-formation 



^ The yield parameter y should not take arbitrary values since it 
is constrained by the IMF and the theory of stellar nucleosynthesis 
(albeit subject to the uncertainties in both). However, it is treated 
in most SAMs as a free parameter. 



3.1 Basic equations of the GCE 

For the purposes of tracing the enrichment of the ISM, we 
still model our galaxies as a single zone with instantaneous 
mixing of gas. We assume that newly produced metals are 
deposited into the cold gas, and may subsequently be ejected 
from the galaxy and mixed with the hot halo gas (or ejected 
from the halo altogether) according to the feedback model 
described above. The metallicity of each new batch of stars 
equals that of the cold gas at the moment of formation. In 
this context, the evolution of the abundance of metals in the 
cold gas is given by 

Gzit) = -mZit)+ezit)+[Gz{t)]^^^^^~[Gzit)]^^^^^^ , (7) 

where G{t) is the total mass of gas and Z{t) is the mass- 
weighted metal abundance, Gz{t) = G{t)Z{t) is the mass 
of gas in the form of metals, i/)(f) is the star formation rate 
and ^{t)Z(t) represents the rate at which metals are de- 
pleted from the ISM by star formation, ez(t) is the rate of 
ejecta of enriched material by dying stars (integrated over 
stellar mass), and the last two terms represent the infall of 
cooled halo gas into the galaxy and the outflow of reheated 
gas from the galaxy. Here we refer generally to 'metals' for 
simplicity, but in fact we apply this equation to each indi- 
vidual element by considering the abundance Zi of a given 
element i instead of the total metallicity Z. For compreheii - 
sive reviews of E q. [T] we direct the reader to iTinslevI (|l9S0l ) 
and IPaeell (|l997h . The modelling of the star formation rate, 
the inflow rate (cooling flows) and the outflow rate (super- 
novae and AGN driven galactic winds) have already been 
sketched in the previous section. The different prescriptions 
shown before relate to the terms in Eq. [7] in the following 
way: 



ip(t)Z{t) = rh^^Zc 



Zh\ 



(mrh + m^gn)Zc 



(8) 



where Zc and Zh are the abundances of the cold ISM gas 
and the hot halo gas, respectively. 

In most SAMs previous to this work, chemical evolu- 
tion was traced in a very simple manner by assuming a 
constant 'effective yield', or mean mass of metals produced 
per mass of stars, and the value of this effective yield was 
treated as a free parameter. In the models presented here, 
we have implemented detailed calculations for the produc- 
tion of heavy elements and the chemical enrichment of the 
ISM (the second term in Eq. [T} in a similar fashion to the 
mo dels oflMatteucci fc~G ibson (19^. lTimmes et all (Il995l ) 
and lPipino fc Matteuccil (j200j). In this framework, not only 
do we trace the evolution of the total metallicity, but we also 
track the distinct elements as well. At this moment we can 
follow the evolution of the abundances of 19 different ele- 
ments, but here we will only discuss a-elements and Fe. By 
Q-elements we mean the composite abundance of N, Na, Ne, 
Mg, Si and S. At any given time, the rate at which an el- 
ement i restored into the interstellar medium is calculated 
according to the following formula 
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ez, = / i'{t-TM)Qrm{t-rM)(t>{M)dM 

I Ml 

+ A 4>{Mb) 

Mb,„ 

^0.5 

f{fl)lp{t - TM2)Qmi{t — rA/2)d/i 



X 



+ (l-A) ip{t-rM)Qm^{t-rM)4>{M)AM 

J Mb^ 

I- Mir 

+ / i^{t-TM)Q^^it-TM)<t>iM)dM (9) 

where Ml and Mu are the lower and upper masses of stars 
formed, 'ip{t) is the star formation rate as before, <I>{M) is the 
initial mass function (IMF), /(/i) is the distribution function 
for the mass of the secondary star in a binary pair (/i = 
M^IMb), tm is the lifetime of a star of mass A/, and Qmi{t) 
represent the fractional mass of an element i restored by a 
star of mass M in the form of both newly synthesised and 
unprocessed material. Although not explicitly dependent on 
time, the quantities Qmi depend on metallicity, which of 
course evolves in time. 

Each of the integrals in Eq. |9] represents the contribu- 
tion to the enrichment by stars in different mass ranges. 
The first integral indicates the contribution of single stars 
with masses between Ml ~ 0.8Mq (the minimum mass 
which can restore gas to the ISM within a Hubble time) and 
MB„^ = 3Mq (the minimum mass of a binary system which 
can give rise to a SNe la event). These stars eject their chem- 
ical by-products through stellar winds and end their lives as 
white dwarfs. The second integral refers to the contribu- 
tion from type la SNe, assuming that these events originate 
from C-0 white dwarfs in binary systems exploding by C- 
deflagration after reaching the Chandrasekhar mass. This 
implies a maximum primary mass of 8Mq, and therefore 



Me 



16Mq and Hn 



Ml 



The pa- 



A/b ' Mb 

rameter A represents the fraction of binary systems with 
total mass in the appropriate range that actually give rise 
to a SNe la event. In essence, A is a free parameter. Chemi- 
cal evolution models of the Milky Way constrain the value of 
this parameter to around A ~ 0.04 — 0.05 by ensuring com- 
patibility with the obs erved present-day ra te of SNe I and 
SNe II in our galaxy l|Francois et al.l [20041 ). However, this 
value results in an unacceptably high abundance of Fe and 
Fe-peak elements in our models. Therefore, we allow this 
parameter to take different values (0.015, 0.02, 0.03, 0.05) 
and constrain it a posteriori by comparison w ith abundance 
ratio s and SNe rates (see also the discussion in lde Plaa et al.l 
I2OO7I ) . The distribution function of the secondary mass frac- 
tion is assumed to follow the law 



(1 + 7)m" 



(10) 



with 7 = 2. A complete description of all the quantities 
involved in the computati on of the SNe I rate can be found i n 
iGreggio fc Renzinil (Il983l ) and lMatteucci fc Greggicl (|l986l ). 
Note that in this scheme, it is the mass of the secondary 
star {M2) that sets the clock for the explosion. This implies 
a specific delay time distribution (DTD) for the explosions 



that may not represent reality (see §|4]2]). The third integral 
represents the mass restored by stars in the mass range 3- 
16 Mq which are either single, or, if binaries, do not produce 
a SN I event. These stars end their lives as white dwarfs 
(M < 8Mq) or as SNe II (M > 8Mq). Finally, the last 
term represents the contribution of short-lived massive stars 
(M > 16 Mq) that explode as SNe II. The fact that we take 
into account the lifetimes of the stars implicitly involves a 
time delay for the different enrichment modes (AGB stellar 
winds, SNe I and SNe II) since the integrands in Eq. |9] are 
by definition zero whenever t < tm- 

Before applying the GCE to the semi-analytic model, 
we tested the chemical evolution code separately. We ran 
simulations using only the chemical evolution algorithm and 
compared the results with simple models with analytic solu- 
tions, i.e., a closed box with either constant and continuous 
star formation or a single initial burst of star formation. 
We also compared our results wi th the output from othe r 
well-tested models, namely that of lFenner fc GibsonI l|2003l ). 
Only after achieving satisfactory agreement in these tests 
did we proceed to implement the chemical evolution into 
the SAM. 



3.2 Ingredients of the GCE 

In the previous section, we introduced several fundamen- 
tal quantities that determine the chemical enrichment of a 
galaxy's cold gas reservoir and described them in a qualita- 
tive manner. In what follows, we will point the reader to the 
different studies that quantify the ingredients of this model, 
and the values adopted for the simulations presented in the 
next section. We discuss the implementation in Appendix 
A. 



3.2.1 Initial Mass Function 

The stellar initial m ass functi o n tha t we use is based on the 
parameterization of IChabrieij (|2003l ): 



_4 g— (log m — log mc)^/2 <T^ 



if m < 1 Mq 
if m > 1 Mq 



(11) 



where in the standard Chabrier IMF, x = 1.3, a = 0.69, 
rric = 0.079 Mq, and the proportionality constants take 
the values A = 0.9098 and B = 0.2539 after normalisation 
in the mass interval 0.1 — 4OAf0. Th is IMF d i ffers s ome- 
what from th e standard power laws of lSalpeteil (Il955l ) and 
iKroupa et al. I lfi993i ) often used in the literature. The reason 
for this choice is consistency with the stellar population syn- 
thesis models that are used to predict magnitudes and colors 
in the simulations, since t hose models use the Chabrier IMF 
l|Bruzual fc CharlotI I200I ) . Note that this expression is for 
the IMF by mass. 

We show below that the results with this IMF were 
not entirely satisfactory. We therefore explored different 
slopes (x = 1.2, 1.1, 1.0) and upper mass limits [Mu = 
60, 100, 120 Mq^. The values of the normalisation constants 
A and B for different values of x and Mjj are computed by 
requiring that 



Mjj 



Ml 



(m) dm : 



(12) 
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We note that these alternate val ues for the slope are within 
the observational uncertainties (|Chabrieij|2003l ), and cer- 
tainly do not represent a radical departure from the observed 
local IMF. 



3.2.2 Stellar lifetimes 

The adopted relation between the evolutionary lifetimes and 
the stellar mass is th at of Padova tracks for solar metallicity 
llBerteUi et al.lll994l '). In principle the lifetimes depend not 
only on mass but also on metallicity. Nevertheless the differ- 
ence in stellar ages for different metal abundances is smaller 
than the age binning in the grid that stores the star forma- 
tion history of the galaxies in our simulations (see Appendix 
A). 



3.2.3 Stellar yields 

Stellar yields are the amount of material that a star can pro- 
duce and eject into the ISM in the form of a given element, 
and are clearly one of the most important ingredients in any 
chemical evolution model. These yields are the quantities 
Qmi in the equations above. In this work we adopt different 
nucleosynthesis prescriptions for stars in the different mass 
ranges. 

Low and intermediate mass stars (0.8 < M/AIq < 8) 
produce He, C, N, and heavy s-process element^ which 
they eject during the formatio n of a planetary nebulae. The 
yields that we adopt are from iKarakas fc Lattanzid l|2007l . 
hereafter KL07). 

Massive stars (M > 8A/q) produce mainly a— elements 
(O, Na, Ne, Mg, Si, S, Ca), some Fe-peak elements, light 
s-process elements and r-process elements. They explode 
as core-collapse Type II SNe. We adopt the yields from 
IWooslev fc WeaveJ (|l995l , hereafter WW95). Note that the 
upper mass limit in this study is 40 A/q . 

Type la SNe are assumed to be C-0 white dwarfs in bi- 
nary systems, exploding by C-deflagration after reaching the 
Chandrasekhar mass via accretion of material from the com- 
panion star. They mainly produ ce Fe and Fe-peak ele ments. 
The yields we adopt are from iNomoto et al.l l| 19971 . here- 
after N97), model W7. When calculating the contribution 
of SNe la, we assume that the primary star also enriches the 
medium as a normal AGB prior to the SN event. 

Except for the SN la yields, which are given only for so- 
lar metallicity, we use metallicity dependent yields, namely 
those tabulated for Z = 0.0002, 0.004, 0.02, interpolat- 
ing when necessary but never extrapolating. Whenever the 
metallicity falls below or above the limiting values, we use 
the yields corresponding to the minimum or maximum Z 
respectively. 

Note that u nlike in some recent st u dies of galactic chem- 
ical evolut i on (iFrancois et al. 2004 : iPipino fc Matteuccil 
|2004 I2OO6I : iNagashima et all l2005al lbl) we do not alter the 
yields in any way in our standard model. We want to see 
if we can fit the data with as few degrees of freedom as 
possible. 



3.3 Delay Time Distribution formulation for SNe 
la 

As mentioned before, and shown in § 14.21 the SN la model 
described in the previous subsection does not seem to be 
the best representation of this phenomenon. In order to test 
other models, we implemented the delay-time -distribution 
(DTD) formalism developed bv lGreggid (|2005| ). In this sce- 
nario, the SN la rate is described by 



Rla{t) = ka 



A{t~T)ip(t~r)DTD{T)dT (13) 



where tl){t) is the star formation rate, r,; is the minimum 
delay time for the SNIa events which we assume to be equal 
to the lifetime of an 8 Mq star, is the maximum delay and 
equal to the lifetime of a 0.8 Mq star, and ka is the number 
of stars per unit mass in a stellar generation defined by 

pMu 

ka = (f>(Tn)m~^ dm. (14) 

■J Ml 

Finally, A{t — r) is the fraction of binary systems which give 
rise to Type la SNe and may, in principle, evolve in time, 
but here we will assume it to be constant. It should be noted 
that in this case, A{t — r) is the fraction relative to the full 
mass range defined by the IMF {Ml — Mu) and not only 
the mass range 3 — 16 Mq as before. To ease the comparison 
with the previous model we will define A{t — r) = Af^-ie, 
where /3-16 is the fraction of stars in the 3 — 16 Mq mass 
range (defined by the IMF) and A has the usual meaning. 

This formulation allows for different SNe la models de- 
pending on the DTD(t) used. In particular w e have cho- 
sen the distribution favoured bvlMan nucci et al] ||2006|). a nd 
parametrized it in similar way as Ma tteucci et al.l ( 20061 ): 



log_DrD(r) = 



1.4-50(logr + 1.3)^ ifT<To 
-0.71 -0.9(logr-f 0.3)^ if r > To 



(15) 



where delay time r is in Gyrs, To = 0.0851Gyr, and the 
distribution function is normalized so that 



DTD{r) dr = 1 



(16) 



^ At this point we do not trace s-process or r-process elements, 
but they will be included in future versions of the code. 



In Figure[T]we show the behaviour of the implemented DTD. 
Note that for a single burst of star formation, about half of 
all SN la explosions occur within the first 100 Myr. 



4 RESULTS 

In this section we present the first results of our model and 
compare them with stellar population studies of a variety of 
early-type galaxies in the local universe. For this purpose, 
we ran simulations for a grid of dark matter haloes of differ- 
ent masses, ranging from 10^^ to IO^^A/q, using both the 
original (instantaneous recycling) and new (full GCE) ver- 
sion of the semi-analytic code. As mentioned before, we limit 
our analysis to the central galaxies of each DM halo. We se- 
lected early-type galaxies from our simulations according to 
the ratio of their bulge-to-total luminosity in the B-band; 
namely, we consider a gala xy to be an early-type when th is 
ratio is larger than 0.4047 l|Simien fc de Vaucouleurslll986l '). 
This selection encompasses both elliptical and SO galaxies. 
Unless otherwise noted, the model results presented in this 
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Figure 1. The normalized DTD proposed by iMannucci et all 
1120061 ) ■ and adopted in this work. 



section are always for central early-type galaxies. We begin 
with a comparison of the results of our new model with those 
of the SOS SAM. We then proceed to compare the predic- 
tions of the new model with observations. For the purpose 
of this comparison, we will only show model galaxies with 
masses above 10^ Mq since the formation history for galax- 
ies below this mass cannot be accurately resolved given the 
mass resolution of the dark matter trees (see §2.1). 



4.1 Impact of the new GCE modelling in galaxy 
observables 

The new chemical evolution modelling affects the physics 
in the SAM in at least three ways: changing the metallic- 
ity of the hot gas changes the cooling rates; changing the 
metallicity of the cold gas changes the amount of dust, and 
therefore observed colours and magnitudes; and the metal- 
licities themselves and their evolution change as well. The 
original version of the SOS SAM did very well at reproducing 
several key properties of galaxies and it is important to ver- 
ify that this is still the case after implementing the detailed 
chemical evolution model. For this purpose we will compare 
simulations with three different 'flavors' for the models: the 
original SAM from SOS, a SAM-I-GCE with standard pa- 
rameters (i.e. a; = 1.3 and A = 0.05), and a SAM-I-GCE 
with parameters that best fit [best-fitting) the abundance 
ratios and metallicities of observed galaxies (x — 1.1 and 
A = 0.015, see §[42}. 

A fundamental feature that any good galaxy formation 
model must reproduce is the luminosity function or stellar 
mass function of galaxies. Given the distribution function 
of dark matter halos and sub-halos predicted by CDM, the 
relationship between stellar mass and halo mass implies a 
specific stellar mass function. The required stellar mass to 
halo mass relationship, in the form of the fraction of baryons 
in the halo that are converted to st ars in a galaxy , has b een 



derived bv lWang etahl l|2006l ) and lMoster et al.i l|2009l ). In 



Figure |5] we show the stellar mass fraction as a function of 
halo mass for the SAM with and without G CE and com- 
pared with the empirical relation obtained bv lMoster et al.l 




11.5 12 



13.5 



Figure 2. The fraction of baryons in the form of stars as a func- 
tion of halo mass for model galajcies in the SAM-I-GCE with the 
best-fitting parameters {x = 1.1 &c A = 0.015, black squares), 
with standard parameters (x = 1.3 ii A = 0.05, green squares), 
and in the original SAM without GCE (red squares). The blue 
lin es mark the empiric al relation, and l-cr uncertainties, derived 
bv lMoster et al.l ||2009|) . 



("2009"). The agreement between all 'fiavors' of the models 
(SOS, standard parameters, best-fitting) and the observa- 
tions is excellent. This implies that the stellar mass function 
in our new models will be nearly identical to that presented 
in SOS, which was shown to agree well with observations. 

In Figure [31 we show the average star formation his- 
tories (SFH) of our model early-type galaxies, for galax- 
ies with different present-day stella r masses. As pointed 
out in previou s studies (see, e.g., iDe Lucia et al.l l2006l : 
SOS: Trager fc Sonicrvillc 200^, hereafter TS09), SAMs with 
AGN feedback do qualitatively reproduce a downsizing- 
like trend (i.e., the higher mass galaxies have shorter SF 
timescales and less extended SFH). This is still true in our 
new models. 

Another very important observational quantity to re- 
produce is the colour-magnitude diagram (CMD) of galax- 
ies. In Figure [4] we show magnitudes and colours of early- 
type galaxies in the SAM with and without GCE. Here there 
is one caveat regarding the calculation of the galaxy lumi- 
nosities. The stellar population models that we use to pre - 
dict the colours and mag nitudes l|Bruzual fc CharlotI [ioosh 
use a fixed standard Chabrier IMF, while we allow the slope 
of this IMF to change when calculating the chemical evo- 
lution. Although this is not self-consistent, such a minor 
change in the IMF parameters should not significantly af- 
fect the predicted colours or magnitudes since the early- 
type galaxies studied here are dominated by old popula- 
tions for which high-mass stars are of little importance (C. 
Conroy, private communication). We divide the CMD into 
red and blue region s using the magnitude-dependent cut of 
iBaldrv et al.l (120041 ). The galaxies form two clear groups: the 
majority in a bright red sequence and a few in a fainter blue 
cloud. The "original" and "best fit" models agree quite well 
with the observed CMD, while the luminous galaxies in the 
"standard model" are slightly too blue. 

In the previous version of the SAM, the fraction of cold 
gas relative to stars in galactic discs at the present time 
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Figure 3. Smoothed average star-formation histories for model 
early-type galaxies in (from bottom to top) the SAM+GCE with 
the best-fitting parameters (x = 1.1 & A = 0.015), with standard 
parameters {x = 1.3 & A = 0.05), and in the original SAM 
without GCE, binned by the stellar mass of the galaxy at 2 = 0. 
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Figure 4. Colour-magnitude diagram (CMD) for model early- 
typo galaxies in the SAM+GCE with the best-fitting parameters 
(x = 1.1 & j4 = 0.015, black stars), with standard parameters 
{x = 1.3 &, A = 0.05, green stars), and in the original SAM with- 
out GCE (red stars). The blue lines mark, from top to bottom, 
the l ocus of the red seque nce, the green valley and the blue cloud 
from lBaldrv et al.l l l2004i) . 
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Figure 5. The gas fraction for model galaxies in the SAM-f GCE 
with the best fitting parameters (x = 1.1 & A = 0.015, black 
symbols), with standard parameters [x = 1.3 &z A = 0.05, green 
symbols), and in the original SAM without GCE (red symbols) 
as a function of u — r colour. Stars and triangles depict early- 
type and disc galaxies, respect ively. The thick b lue line marks 
the median of the sample from iKannappanI ||2004| ). and the thin 
lines mark the I — a deviation. 



was used to calibrate the models by comp arison with the 
observational estimates of lBell et al.l l|2003l ) for morphologi- 
cally late- type galaxies (see Figure 5 in SOS). This property 
is well reproduced by all test cases after implementing the 
GCE. Here we also compare the gas fraction o f our galax- 
ies, bo th early-type and discs, with those from IKannappanI 
l)2004h as a function oi u — r colour, as shown in Figure [5] 
The agreement in the slope, scatter and zero-point of the re- 
lation is quite good for all models, especially when the new 
model of galactic chemical evolution is included. 

The first evidence for the mass-metallicity relation can 
be seen when looking at the metallicity distribution function 
of galaxies of difi'erent masses, which we show in Figure[S]for 
our three test cases. All of the models agree qualitatively, 
showing an increasing mean of the distribution as the mass 
range increases. At a given mass, however, the distributions 
for the best-fitting model are shifted to higher metallicities 
since galaxies in this simulation are more metal rich (see 
below) . 

Summarising, we have seen that including a detailed 
chemical evolution model in the SAM has a minor effect on 
the predicted formation histories and present-day properties 
of galaxies, and therefore does not require a re-calibration of 
the free parameters of the model. In the next subsection we 
will investigate the predicted metallicities and abundances, 
which are indeed affected by the GCE. 



4.2 Chemistry of Early-type Galaxies 

The main effect of the new treatment of chemical evolu- 
tion is reflected in the metallicity and abundance ratios of 
the galaxies. Most SAMs reproduce fairly well the mass- 
metallicity relation of galaxies (with effective yield treated 
as a free parameter), but to date they have been unsuc- 
cessful in fitting the slope of the mass-[a/Fe] relation (e.g.. 
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Figure 6. Average metallicity distributions for model early-type 
galaxies in (from bottom to top) the SAM+GCE with the best 
fitting parameters {x = 1.1 h A = 0.015), with standard parame- 
ters (a; = 1.3 & A = 0.05), and in the SAM without GCE, binned 
by stellar mass. 



iNagashima et all l2005bl : IPipino et al] |2008| ). This is the 
main challenge that we address in this study. 

The galaxy sample used for comparison is that de- 



scribed in iTrager et all (|2000al ). This sample has been re- 
analysed using the updated stellar pop ulation synthesis 
method presented in iTrager et al.l l|2008h . which is sensi- 
tive to age, metallicity and abun dance ratios. In the current 
study, we use models based on the lBruzual fc Chariot! (|2003l ) 
mod els with index var iations due to abundance ratios taken 
from iLee et al.l l|2009l ). Inferred stellar population parame- 
ters are tabulated in Appendix B. When making the com- 
parisons, we use the stellar mass of the simulated galaxies 
and the inferred dynamical mass of the observed ones. This 
should not introduce a significant bias since the dynamical 
mass is a good tracer of the stellar mass within one effectiv e 
radius for most early-type galaxies ( CappcUari et al1l2006l ). 
The abunda nces presented here w ere normalised to the solar 
values from iGrevesse et "all (Il996l ). We have also computed 
the present epoch SNe la and II ra tes for our galax i es and 
com pared them w i th the results of iMannucci et al] (|2005l ) 



com pared them w i th the 
and lSuUivan et all l|2006h 



4-2.1 Dependence on IMF slope and SN la fraction 

In Figures[7|and[8]we show the relation between total metal- 
licity ([Z/H]) and stellar mass (M«). From these figures, we 
see that in the "old" SAMs, galaxies tended to be too metal- 
poor compared to the observations (as seen in TS09), and 
implementing the detailed chemical evolution with the stan- 
dard IMF improves the results only slightly. We explore the 
parameter space of the GCE equations to see if it is possible 
to improve the agreement with the observations. Specifically, 
the parameters allowed to vary are the fraction of binaries 
that give rise to SNe la (A in Eq. ^ , the slope of the IMF 
above 1 Mq {x in Eq. Ilip and the upper mass limit of the 
IMF (Mc7 in Eq.[5|. However, we will only show the results 
for Mu = 40 Mq for the following reason. Given that our 
chosen SN II yields (WW95) are only tabulated up to that 
value, we are forced to assume that the yields relative to 
the initial mass remain constant and equal to those of a 
40 Mq star for stars above this value when running simu- 
lations with higher upper mass limits. This proves to be an 
unreliable assumption, as reflected for example in the non- 
monotonic behaviour of the [a/Fe] ratio as a function of the 
stellar mass of our galaxies, in serious disagreement with 
observations. Other properties, such as the SNe la rate, are 
not significantly affected by this parameter, as expected. 

We also note that, even though we have varied the slope 
of the IMF, the supernova feedback efficiency remains the 
same because the prescription for the SN feedback energet- 
ics and the SN chemical enrichment are decoupled; the SN 
feedback efficiency is set manually as a free parameter (see 
Eq. [3]), independent of the IMF. In this sense, the model 
is not entirely self-consistent. However, this choice makes 
it easier to interpret the effect of changing the IMF in the 
models. 

We now compare the predictions for the observed galaxy 
sample, the original SAM, and the new SAM with GCE 
with a standard Chabrier IMF, a shallower IMF (a; = 1.1) 
and two choices of the parameter A (0.015 and 0.05). The 
use of a shallower IMF results in an upwards shift in the 
zero-point of the mass-metallicity relation, bringing it into 
much better agreement with the observations. This increase 
in metallicity with a shallower IMF is expected since more 
massive stars are produced and therefore the gas is enriched 
more efficiently. The metallicities of the model galaxies show 
almost no dependence on the parameter A, which is not 
surprising since this parameter mainly controls the ratio of 
type la to type II supernovae, affecting abundance ratios but 
not the overall metallicity (in essence the production of a- 
elements and Fe-peak elements compensate for one another). 

A primary advantage of our new model is that we can 
now calculate abundance ratios for our galaxies, which gives 
us yet another property to compare with observations and 
set further constraints on the models. This is shown in Fig- 
ures |5] and [TOl where we plot the [a/Fe] ratio against stellar 
mass for the same parameter choices as before. For the sim- 
ulated galaxies, we consider the abundance of a-elements to 
be the com posite abundance of N, Na, Ne, Mg, Si, and S (cf. 
ITrager eTa l. 2000a) . If we assume a value of the parameter 
A ~ 0.05 commonly used in the literature, the abundance ra- 
tios are far too low. The overall values of [a/Fe] can be raised 
by decreasing the value of A (i.e. producing fewer SNe la and 
consequently less iron). However, even with this decreased 
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Figure 7. The relationship between metallicity and stellar mass 
for the galaxies in our simulations and in the observational 
sample of TOO. Symbols — pink squares: original SAM; red 
crosses: SAM+GCE with x = 1.3 and A = 0.015; blue crosses: 
SAM+GCE with x = 1.3 and A = 0.0 5; black stars with error 
bars — galaxies from lTrager et al.l ll2000al) . reanalysed as described 
in the text. Note the poor agreement of the model galaxies with 
the observations in all cases. 
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Figure 9. Relation between the [a/Fe] ratio and stellar mass. 
Symbols - red crosses: SAM+GCE with x = 1.3 and A = 0.015; 
blue crosses: SAM+GCE with x = 1.3 and A = 0.05 ; black stars 
with error bars: galaxies from iTrager et al. I l l2000al'l . reanalysed 
as discussed in the text. Note again the poor agreement of the 
model galaxies with the observations. 
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Figure 8. As in Figure [J] except with a shallower high-mass 
slope for the IMF {x = 1.1) used in the SAM+GCE. Note the 
significantly better agreement of the model galaxies with the ob- 
servations. 
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Figure 10. As in Figure[9l except x = 1.1 for the model galaxies 
from the SAM+GCE. Note the better agreement with observa- 
tions when using a low value of the SN la fraction A. 



A and the standard IMF, we see that the relation is too flat 
or even has a slightly negative slope, and for the galaxies at 
the high mass end the [a/Fe] ratio is insufflciently enhanced. 
On the other hand, a shallower IMF, combined with a lower 
value of A, produces model galaxies in far better agreement 
with observations. The flatter IMF increases the slope of the 
relation while a low A brings up the zero-point. Moreover, it 
is interesting to note that lowering A also increases the slope 
(at any flxed IMF). This small steepening of the [a/Fe]-mass 
relation is due to the metallicity dependence of the yields, 
since the higher the initial metallicity, the higher the [a/Fe] 
in the yields. A lower fraction of SNe la (and consequently 
more type II) implies a slightly faster overall enrichment, 
and therefore galaxies spend more time forming stars in a 
regime of enhanced [a/Fe] (higher metallicity yields). 



In summary, we had expected that the inclusion of AGN 
feedback in the semi-analytic models might solve the prob- 
lems that previous studies have encountered in trying to 
reproduce the observed trend between mass and [a/Fe] ra- 
tio, because the quenching due to AGN does lead to more 
massive galaxies having shorter formation timescales in the 
models. However, apparently the effect of "downsizing" on 
the trend of [a/Fe] is very small and a flatter IMF is re- 
quired to achieve agreement between the models and the 
observations. This does not undermine the potential impor- 
tance of AGN feedback, which appears to be a promising 
mechanism for solving many of the other problems experi- 
enced by earlier generations of models (such as the over- 
cooling problem). In any case, it is encouraging that with 
minor variations (within the observational uncertainties) in 
the chemical evolution parameters, we can for the first time 
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Figure 11. Present day SN la rate as a function of specific 
SFR. Here we show all model galaxies, regardless of morphol- 
ogy. Red stars are SAM+GCE with x = 1.3 and A = 0.015; blue 
stars SAM+GCE with x = 1.3 and A = 0.05; black squares are 
observati ons from ISuUivan et ah 1 1 I2OO6I) and black crosses from 
[Mannuc ci et al.l ||2005|) . The conversion of galaxy type into spe- 
cific SFR for the M annucci et al. data points is the same as in 
ISulhvan et al] l l200d) . 



obtain very good agreement with the observed mass-[Q:/Fe] 
in a semi-analytic model. 

One concern is that the comparison between our mod- 
els and the observations is not strictly rigorous since we 
are showing stellar mass-weighted abundances for the mod- 
els, while stellar population studies derive abundances from 
line-strength indexes from integrated spectra, which are 
themselves light-weighted quantities. However, TS09 have 
shown that the SSP-equivalent (absorption-line-weighted) 
metallicity correlates very well with its mass-weighted and 
light-weighted counterparts. In future work, nevertheless, 
we will synthesise line strengths for the galaxies in our 
simulations and calculate abundance ratios in the same 
way as is done in the observational data. One additional 
worry is the effect of limited aperture size on the compar- 
ison, as early- type galaxies are well-known to have signif- 
icant line-strength gradients. These gradients imply how- 
ever mild metallictty gra dients but no abundance ratio gra- 
Davies et all 1 19931 : iMehlert et all 



dients whatsoever (e.g.. . . , . 

l2003l : ISanchez-Blazauez et al.ll2007l .TS09'). Therefore we are 
confident that trends in [a/Fe] with mass are trustworthy. 
We note here (as described in Appendix B) that the data 
plotted in the figures have been constructed to appear as 
if the galaxies were at the distance of the Coma cluster 
and observed thro ugh fibre apertures of diameter 2'.' 7 (see 
iTrager et al.|[20o3 ). While this does not eliminate gradient 
eflects on the inferred metallicities, it reduces their magni- 
tude to an ofltset of roughly -0.1 dex (TS09). 

Supernovae rates provide further independent con- 
straints our models. We calculate the predicted supernovae 
rates for ou r model galaxie s and compare them with those 
derived by Mannucci et all (j2005) and Su llivan et al.1 (120061 ) 
for a large sample of galaxies in the nearby universe. In Fig- 
ures [TT] and [T5] we show the present day type la supernova 
rate in units of SN events per year per unit stellar mass ver- 
sus the specific star formation rate (SSFR, star formation 
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Figure 12. As in Figure 1111 except using the shallower IMF 
(x = 1.1) for the model galaxies from the SAM-I-GCE. 
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Figure 13. Present-day SN II rate as a function of SSFR. Sym- 
bols are as in Figure [TT] Results are shown for the standard IMF 
{x = 1.3). We remind the reader that the upper-mass limit on the 
IMF is 40 Mq 
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Figure 14. The same as Figure [131 except we use the shallower 
IMF slope (a: = 1.1) for the model galaxies from the SAM-I-GCE. 
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rate per unit stellar mass), and in Figures [13] and [14] we show 
the rates for Type II SNe. Here we show all the model galax- 
ies regardless of morphology, since the early-type galaxies 
populate only the lower SSFR side of these diagrams, and 
both comparison samples include galaxies of all morpholog- 
ical types. The galaxies in these sample are from a mixture 
of field and small cluster environments, but so are our model 
galaxies. 

The slope of the IMF has little effect on the predicted 
SNIa rates; only the A parameter has a significant infiuence 
on the results. From Figures [TTl and fT^ we see that no combi- 
nation of IMF slope and fraction of binaries that yield SNe 
la can fit the observations over the whole range of SSFR. 
However, it is interesting to notice that SNe la rates of star 
forming galaxies (SSFR > 10"^"'^) are very well matched 
by models with a high value of A, while a low value of A is 
a better match for passive galaxies (SSFR < lO"^"'^). This 
behaviour, which is seen regardless of the slope of the IMF, 
is alm ost certainly due to the ch osen supernova la model. 
In the iGreggio fc Renzinil (|l983D formalism which we use, 
the Delay Time Distribution (DTD) of the SN la explo- 
sions is given by a convolution of the distribution of sec- 
ondary masses (in binary systems) and the lifetime of the 
secondary sta r. On the other hand, from the same obser- 
vational data, iMannucci et al.l (|2006l ) derived a DTD with 
two components: a prompt peak and a later plateau, each 
encompassing half of the SNe la. Other authors have also 
reached similar conclusions about the delay-time distribu- 
tion of type la explo sions (e.g. lScannapieco fc Bildstenll2005l : 
iDahlen et alll 20081 ). This bimodal DTD effectively enhances 
the production of type la supernovae in star forming galax- 
ies, exactly where a higher fraction of SNIa are needed in 
our models. We therefore expect that using a two-population 
DTD with a significant prompt component will alleviate the 
differences between our model predictions and the observed 
Type la rates, and we show that this is the case below. 

The Type II SN rates, on the other hand, show very 
good agreement with the observations over the whole range 
of SSFR. In this case, all variations of x and A give the same 
qualitative results, although the models with the "standard" 
parameter values are slightly better. Nevertheless, as we 
have shown above, not all combinations of IMF slopes and 
values of A (SN la producing fraction of binaries) produce 
model galaxies that agree with observed metallicities and 
abundance ratios. It is only for those models with a shal- 
lower IMF [x ~ 1.1) and lower SNe la fraction among bina- 
ries [A ~ 0.01-0.02) that we can reasonably reproduce the 
full set of observations. 

To summarise, after implementing detailed chemical 
evolution in the semi-analytic model of SOS, we can now 
reproduce the mass-metallicity and mass-[a/Fe] relations 
for local early-type galaxies, provided we use a slightly flat- 
ter Chabrier IMF and a low fraction of binaries giving rise 
to SN la. The predicted rates of Type la SNe show a strong 
dependence on the fraction of binaries that yield such an 
event (our parameter A) but not on the slope of the IMF, 
with the rates in star forming galaxies better matched by a 
high value of A, while those in passive galaxies are better 
with a low value of A. However, for a standard IMF and 
high value of A, the galaxies are a bit too metal-poor and, 
most importantly, the abundance ratios are extremely low, 
in severe disagreement with observations. 



4-2.2 Bimodal Delay Time Distribution for Type la 
Supernovae 

In the previous section, we speculated that a bimodal distri- 
bution with a prompt population of SNe la explotions would 
give a better match for the observed supernova rate s. Given 
the a nalytical nature of the DTD formulation (jGreggiol 
l2005h . it is fairly straightforward to implement and test 
thi s hypothesis. We have implemented the DTD proposed 
by iMannucci et al.l (|2006l ). Figure [15] shows the four ob- 
servational constraints used to test the models (metallic- 
ity vs. stellar mass, [a/Fe] vs. stellar mass, type la SNR 
vs. SSFR, and type II SNR vs. SSFR). Clearly the new 
double-peaked SN la DTD model gives a better match to 
the Type la SN rates, while maintaining the good agreement 
in the other galactic properties. However, for this model, 
the best-fltting parameters for the IMF slope and the SNIa 
binary fraction are slightly different; specifically, x = 1.15 
and A = 0.03. This value for the I MF is in even better 
agreement with some rec ent studies (Bald rv fc Glazebrookl 
[2OO3I : IWilkins et al.ll2008l ) then our previous 'best' value of 
X = 1.1. 

Finally, in Figure [TH] we show the abundance ratios of 
some individual elements for the galaxies in our best fit- 
ting models using the classic SNe la recipe (x = 1.1 and 
A = 0.015) and the bimodal DTD {x = 1.15 and A = 0.03). 
With the exception of C and N, which are slightly higher for 
the latter, all the elements follow the same trends in both 
models. In particular Mg, even though it is under-abundant 
with respect to other a-elements, is the one element that 
best follows the observed trend of [a/Fe]. This leads us to 
believe that the abundances derived from the stellar popula- 
tion analysis may be predominantly driven by Mg, and may 
not necessarily reflect of all the a-elements. There is also an 
excess of Ni in the Fe-peak group and a decreasing trend of 
[C/Fe] wi th increasing galactic mass wh i ch is apparently not 
observed (jSanchez-Blazquez et al.|[2003l : [Graves fc Schiavonl 
|2008l). This should not be considered a flaw of the model, 
since the abundances of individual elements are very sensi- 
tive to the chosen yields. We expect that future line-strength 
observations, interpret ed with next-generation st ellar popu- 
lation models (such as lSchiavonll2007l : iLee et al. l 2009), will 
provide an interesting test of our models, including the set 
of assumed yields (KL07 + WW95 + N97). 



4-2.3 Dependence on other model parameters 

We have also explored whether we could match the obser- 
vations investigated above by varying the galaxy formation 
parameters of the SAM instead of the IMF and binary frac- 
tion parameters. For this purpose, we ran several simulations 
in which we modifled the star formation efficiency {Ak in 
Eq. [2]|, the SN feedback efficiency (eo^ in Eq. [3} and the 
virial velocity below which ejection of reheated gas from the 
halo into the diffuse intergalactic medium becomes impor- 
tant (Vcjcct in Eq.[3]). When the star formation efficiency was 
increased by a constant factor, a few more massive galax- 
ies were produced and the slope of the mass- [a/Fe] rela- 
tion increased slightly but the slope and zero point of the 
mass-metallicity relation did not change. A higher SF effi- 
ciency implies that the cold gas is consumed more rapidly 
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Figure 15. Clockwise, starting from the top left panel: (A) [Z/H] vs. stellar mass; (B) [a/Fe] vs. stellar mass; (C) Type la SNR vs. SSFR; 
(D) Type II SNR vs. SSFR. Symbols - red crosses: SAM+G CE with bimodal DTD for SN e la {x = 1.15, A = 0.0 3); T OO: reanalysed 
metallicities and abundance ratios from iTrager et al] (|20003), M05 and S06: SN rates from iMannucci et al.l | |2005^ and ISuUi van et al.l 
||2006| ') respectively. 



and therefore the timescale for star formation is shorter, 
which is why the mass-[a/Fe] relation was affected. 

If Ak was allowed to increase with increasing galac- 
tic baryonic mass, a considerable number of very massive 
galaxies were produced and the slope of the mass-metallicity 
relation increased mildly (but not the zero point). On the 
other hand, if the factor decreased with increasing galactic 
mass, the trends remained the same but no galaxies above 
10" Mq were produced. This excess or lack of high mass 
galaxies arises because star formation in the biggest systems 
is either boosted or suppressed by this mass-dependent vari- 
ation of the SF efSciency. 

The effects of reducing the SN feedback efficiency and 
the T4joct parameter were roughly the same. In both cases 
the slope and zero point of the mass-metallicity relation 
increased slightly and the central galaxies of the DM haloes 
were on average more massive, but the mass-[Q/Fe] relation 
did not change and an unrealistically large number of low- 
mass satellite galaxies was also produced. This excess of low- 
mass satellites is due to the fact that small galaxies retain 
their gas more efficiently when these parameters that control 
the SN feedback are decreased. 

Overall, the effect of changing these other parameters 
is not as strong as flattening the IMF, and also destroys the 
agreement with other well-calibrated observations, such as 



the luminosity function, the metallicity distribution func- 
tion, and the cold gas fraction. We therefore conclude that 
our results are robust to the values of these free parameters. 

As a final test, we modified some of the yields. Specif- 
ically, we decreased the Fe yield of SNe II by half and in- 
creased the Mg yield by a factor of four. Reducing the Fe 
yield had very little effect, indicating that the bulk of the 
Fe comes from SNe la, as expected. Changing the Mg yield 
raised the zero point of the relations, as ex pected. How- 
ever, n either of these changes affected the slope. [Pipino et ahl 
(|2008l ) reached a similar conclusion about the yields and 
other parameters when exploring the parameter space in 
their models. 



5 DISCUSSION AND CONCLUSIONS 

We have implemented detailed galactic chemical evolution in 
a semi-analytic model, and use the resulting model to study 
the metal enrichment of early type galaxie s in the local uni- 
verse. The base SAM is that presented in ISomerville et al.] 
(|2008h . We take into account the effects of galaxy mergers, 
inflow of cold gas, and SN and AGN driven outflows, as well 
as the production of metals by SN la, SN II and AGB stars. 
Unlike most previous SAMs we discard the instantaneous 
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Figure 16. [el/Fe] versus stellar mass for so me ct-elements and Fe-p eak elements. Only predictions f rom our best- fitting models are 
shown. Symbols - black stars; Classic SN la llGreegio &: Renzinilll98a) : red crosses: DTD formulation llGregeidl2005t) . The parameters 
are {x = 1.1, A = 0.015) and {x = 1.15, A = 0.03) respectively. 



recycling approximation by properly accounting for the fi- 
nite lifetimes of the stars and also make use of metallicity 
dependent yields. 

We run our SAM-fGCE simulations in a grid of dark 
matter haloes ranging over present-day masses of 10^^ to 
10" Mq. We allow the slope of the IMF and the fraction 
of binaries that produce a SN la event to vary, and com- 
pare our results with the observed trends of metallicity and 
abundance ratio ([a/Fe]) against stellar mass of the galax- 
ies, as well as the supernova rate (both type la and II) as 
a function of specific star formation rate. Only the mod- 
els with a shallow IMF {x — 1.1) and a low fraction of SN 
la from binaries {A ~ 0.015) match all four observations 
of early-type galaxies simultaneously. A slightly flatter than 
standard IMF is necessary in order to produce more massive 
stars, which enrich the interstellar medium more efficiently, 
making the galaxies in our simulations become more metal 
rich and improving the agreement with the data. The pro- 
duction of more massive stars, along with the fact that the 
star formation histories are more extended in time as the 
galaxy mass decreases, helps to achieve the correct trend of 
increasing [a/Fe] with increasing galaxy stellar mass. How- 
ever, it is also necessary to invoke a low fraction of SNe 
la to raise the zero-point of this relation. We also predict 
abundance patterns for a variety of elements for early-type 



galaxies at z = in our fiducial model. These predictions 
will be interesting to compare with future observations. 

From studying the SNe la rates, we find evidence sup- 
porting a 'two-population' distribution for the type la explo- 
sions, since galaxies with high specific star formation rates 
are better matched by models with a high fraction of bina- 
ries that explode as SN la (A) while those with low SSFR 
require a low value of A. We tested whether the use of a 
more realistic (bimodal) delay-time distribution of type la 
supernovae would, in fact, improve the results. After im- 
plementing the DTD formulation for SNe la and using a 
bimodal distribution with a prompt peak and an extended 
plateau, we found very good agreement with the SN rates, 
while still matching the trends of [Z/H] and [a/Fe] with stel- 
lar mass, although the best values for the slope of the IMF 
changed slightly and the fraction of SNela binaries needed to 
be doubled. Our favored model is now one with a Chabrier- 
like IMF with a slope of x — 1.15, a SNe la binary fraction 
of ^ = 0.03 (relative to the 3-16 Mq range, A ~ 0.0014 
relative to the full range of masses defining the IMF), and 
a bimodal d elay-time-distrib u tion f or Type la SN events as 
proposed bv lMannucci et ahl (|2006| ). 

We have also studied the effects of varying the galaxy 
formation parameters in the SAM, but found that we were 
unable to reproduce the observations in this way. We there- 
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fore conclude that our results are robust to the values of the 
free parameters in the SAM. 

This is not the fir st time that a GCE mod el has been 
applied within a S AM . iNaeashima et alJ ()2005al lbh have also 
constructed such a model, and their 'superwind' model re- 
sembles our model. They first obtained fairly good agree- 
ment with observations of ICM abundances in galaxy clus- 
ters, however the same models failed to reproduced the trend 
of increasing [a/Fe] with increasing galactic mass. One of 
the main differences between their models and ours is in 
fact the IMF. They use a Kennicutt IMF {x = 1.5) for qui- 
escent star formation and a flat IMF (x — 0.0) for stars 
formed in bursts. This flat IMF is rather extreme, while the 
proposed modification in our models is small, in fact within 
the observational uncertainties. Namely, we require the same 
'shallow' IMF {x = 1.15) for all modes of star formation. An- 
other model that reproduces the obser ved scaling of abun- 
dance ratio with galaxy mass is that of IPipino fc Matteuccil 
l|2004l . [ 2006h . However they consider a very different sce- 
nario for galaxy formation, the monolithic collapse scenario, 
and allow for galaxy mergers only in the form of asec- 
ond infalling episode. More recently, IPipino et al.l (|2008l l 
also coupled GCE to a SAM, but they also failed to match 
the mass-[a/Fe] relation. They claim that flattening the 
IMF can not solve this problem, in contradiction with our 
findings. It is worth mentioning that none of these mod- 
els include AGN fe edback; only the GallCS model used by 
IPipino et al.l l|2008h has some form of halo quenching sim- 
ply by shutting down the flow of cold gas onto galaxies with 
masses larger than 10^^ Mq. However, contrary to our ex- 
pectations, we find that SF quenching by AGN is not a key 
factor in our success at reproducing the mass-[Q!/Fe] rela- 
tion, even though the AGN feedback in our models leads 
to shorter formation times for the more massive galaxies. 
We find that a slight flattening of the IMF is essential to 
achieve agreement between the model and the observations. 
AGN feedback, nonetheless, is likely to play an important 
role in reproducing other galaxy observations, such as the 
stellar mass or luminosity function and color bimodality. 

Our best-fitting IMF, nonetheless, is consistent with 
observati ons; the slope i s within the observational uncer- 
tainty of IChabried (l2003h and agree s remarkably well with 
the results of Baldrv fc Glazebrookl (jjOOS), who found that 
ultraviolet to near-infrared galaxy luminosity densities re- 
quire an IMF with a slope of 1 .15 =fc 0.2. The same slope was 
found bv lWilkins et al.] (|2008l ) when trying to reconcile the 
redshift evolution of the observed stellar mass density with 
the cosmic SFH using a constant and universal IMF. The 
agreement, however, holds only at low redshift. In a forth- 
coming paper, Wilkins et al. (in preparation) propose an 
evolving IMF as a plausible solution. Such an IMF should be 
strongly top-hea vy at high redshift (Hopkins, private com- 
munication) . . van DokkumI \200^ also suggests an evolving 
Chabrier-like IMF based on comparing the evolution of the 
M/L ratios of early-type galaxies to their colour evolution, 
but in this case the change is in the characteristic mass (nic 
in Eq. Ilip rather than the slope, making the IMF "bottom- 
light" at high redshift. Such evolving IMFs could, in princi- 
ple, work in favour of the trends of [a/Fe] with stellar mass 
and SNR with SSFR since they produce either more SN II 
or fewer SN la progenitors at earlier times when massive 
early-type galaxies create most of their stars. This scenario 



remains to be tested, and moreover the issue of an evolving 
IMF is open to considerable debate given th e large uncer- 
taintics on its constraints. On a different note, iMeurer et al.l 
(2009) has claimed evidence for an IMF that depends on 
galactic surface brightness (or surface density) as a plausible 
explanation for an observed variation in Hc/FUV flux ratio. 
However, they invoke variations that are an order of magni- 
tude larger than the deviation of our best-fittin g slope from 
the st andard value. Finally, very recently, Calura fc Mencil 
(j2009h have claimed that a constant IMF cannot account 
for the trends of [Z/H] and [a/Fe] with velocity dispersion 
in elliptical galaxies and have proposed an IMF with a slope 
depending on the SFR {x — 1.35 for low star-forming sys- 
tems and X = 1 for high star-forming systems) in order to 
explain them. However, their chemical evolution is not cou- 
pled to a semi-analytic model, but is com puted a posterior i 
with SFHs extracted from the SAM of M enci et all (|2008l ). 
and therefore the flows of enriched gas from the galaxies into 
the halos and back again are not tracked, unlike in the our 
model. 

In future work we will apply this model to other ques- 
tions such as the abundances of different components of spi- 
ral galaxies like the Milky Way (e.g. the disk, bulge, and stel- 
lar halo), abundances in clusters vs. the field, abundances 
in the intra-cluster gas, and to the evolution of metals over 
cosmic time. 
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Table Al. Structure of the grid where the star formation and 
abundance histories are stored. All quantities are in Gyr. 
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1.28 



The initial size of the bins is 0.01 Gyr because it is the maximum 
possible value of the time-step in our simulations. 
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APPENDIX A: IMPLEMENTATION 
ALGORITHM OF THE GCE MODELLING 

To compute the chemical enrichment consistently as de- 
scribed in Section 2, it is imperative to keep track of the star 
formation history of each galaxies as well as the metallicity 
of the ISM as a function of time in these galaxies. Given 
the large number of galaxies in the simulations, it is very 
expensive, both in computing time and in physical memory, 
to store these quantities in a linear time grid with reasonable 
resolution. In order to overcome these limitations, we use an 
age grid where first bin of this grid represents t = tnow at 
any moment during the calculations and the size (in time) 
of the bins increase as we go to older ages. The scaling of 
the bins is related to the lifetimes of the stars, so that it 
increases as the lifetimes of the stars become progressively 
larger. The size of the bins only increases, but each subse- 
quent bin is not necessarily larger than the previous one. 
In Table lAll we show the structure of this binning grid. In 
this grid we store the amount of cold gas mass turned into 



stars, AM*, and the metallicity of the cold gas, As the 
systems evolve this information is 'pushed down' to older 
bins, directly if the subsequent bin is of the same size or 
added if the bin is of a larger size. In the later case, the cold 
gas masses are added directly whereas the the metallicity is 
averaged, weighted by the size of the bins. In essence, what 
we extract from this grid are time-averaged quantities. 

It is also important how we handle this information 
when galaxies merge, that is, how we combine the star for- 
mation and metallicity histories of the merging galaxies. 
This is done on a bin-by-bin basis. At a given age bin, AM* 
is added directly and the metallicities are averaged weighted 
by the corresponding AAf*. In mathematical form. 



AM4i\ = AM^\i]+ AM^\i], 

* AM,^[i] + ZfAM^[i 



AM^[{\ + AMP[i 



(Al) 
(A2) 



where A and B are the merging galaxies and i the age bin. 
Summarising, the star formation rate and the cold gas metal- 
licity used in Eq. |9] to calculate the enrichment of the ISM 
are average quantities. 



APPENDIX B: GALAXY SAMPLE 

The sample of galaxies u sed to test our models was taken 
from lTrager et all ll2000ah and r eanalysed with the method 
described in Trager et al. l|2008l '). Specifically, we used the 
iBruzual &: CharlotI ( 200S ) stellar population models and 
modified the line strengths when [a/Fe] 7^ using 'response 
functions' that have been calculated for each star in each 
isoch rone, a significant improvement over pre vious methods 
(e.g., iTripicco fc Beilll995l:l Korn et al ] 120051). More details 
can be found in iTrager et all (120081 ) and iLee et ahl (|2009t ). 
In Table IBll we present the results from this new stellar 
population synthesis analysis. 

We also present dynamical mass computed using 



Mdyn = 465(7 V,, (Bl) 

where a is the central velocity dispersion in kms^^ (here 
we use the velocity dispersion taken through an aperture 
of diameter re/8) a nd re is effective radius in parsecs (see 
iTrager et al]|2000bl . for more details). 



This paper has been typeset from a linjX/ I^I^TjX file prepared 
by the author. 



^ The SFR at the time of interest is computed dividing AAf* by 
the size of the bin. 



GCEin 



Table Bl. Velocity dispersions, ages, metallicities, enhajicement ratios and dynamical masses of the 
sample galaxies. 



Galaxy Type Log(<7) log(Agc) [Z/H] [E/Fc] log(Mdyn) 



ESO358-G06 


SO 


1 


.763 


it 


0, 


.072 


0, 


,807 


+0. 
-0, 


.114 
,250 


— 


.01 +0 

•431+0 


.083 
.098 


— 0, 


.005 


+ (3 
-0 


,066 
,045 


9, 


.24 


it 





.43 


ES0358-G25 


SO 


1 


.690 


± 


0, 


,072 


0, 


,779 


+0. 
-0, 


.326 
,295 


— 0, 


.7321°' 


,159 
,144 


— 0, 


,044 


+0, 
-0 


,126 
,096 


9 


.18 


± 





.39 


ESO358-G50 


so 


1 


.732 




0, 


.072 


0, 


,707 


+0. 
-0, 


.129 
,235 


— 0, 


.110+0, 

.4123g 


,098 
.098 


— 0, 


,076 


+0, 
-0 


,076 
,056 


9, 


.25 


it 


0, 


.48 


ES0358-G59 


so 


1 


.653 


± 


0, 


.072 


0, 


,544 


-f 0. 
-0, 


.235 
,068 


-0, 


•266tg 


.053 
.114 


0, 


,045 


+0 
-0 


,045 
.035 


9, 


.00 


it 


0, 


.45 


ESO359-G02 


so 


1 


.763 


± 


0, 


.072 


0, 


,603 


+0. 
-0, 


.159 
,189 


-0, 


no .1+0 

•904Ig 


.220 
.189 


-0, 


,164 


+0 
-0 


.187 
.500 


9, 


.36 


it 





.46 


1C19d3 


so 


1 


.763 


± 


0, 


,072 


0, 


,397 


+0. 
-0 


114 

,068 


0, 




,083 
,098 


0, 


,085 


+0, 

-0 


.045 
.035 


9, 


.22 


it 


0, 


.50 


IC2006 


E 


2, 


.134 


± 


0, 


,024 


1, 


,126 


+[) 
-0 


()98 
,114 


0, 


1O + 

.049Iq 


068 
,083 


0, 


162 


+(3 

-0 


,025 
,(315 


10 


.30 


± 


0, 


.33 


NGC0221 


E 


1 


.731 


± 


0. 


,021 


0, 


,975 


+0. 
-0, 


.068 
,053 


— 


.398^0 


,038 
.038 


— 0, 


,026 


+0, 
-0, 


,025 
.006 


8, 


.32 


it 


0, 


.28 


NGC0224 


E 


2, 


.185 


± 


0, 


,008 


0, 


,926 


+0, 

-0, 


,129 
,159 


0, 


oni +0 

.2817;q 


.068 
.038 


0, 


,153 


+0 
-0 


.025 
.006 


10 


.61 


it 





.15 


NGC0315 


E 


2, 


.486 


± 


0, 


,005 


0, 


,860 


+0. 
-0, 


.098 
,083 


0, 


^+0 

.316_JJ 


,038 
,023 


0, 


,239 


+0 
-0 


.015 
.005 


11 


.85 


± 


0, 


.35 


NGC0507 


E 


2 


.445 


± 


0, 


,009 


0, 


,601 


+0. 
-0, 


.129 
,098 


0, 


nro + O 

.359_JJ 


(368 
,053 


0, 


179 


+0, 
-0 


.035 
.015 


11 


.89 


± 


0, 


.34 


NGC0584 


E 


2 


.278 


± 


0, 


,005 


0, 


,531 


+0 
-0 


053 
.023 


0, 


01 0+0 
.318+0 


023 
,023 


0, 


172 


+0 

-0 


,015 
,005 


10 


.69 


± 





.30 


NGC0636 


E 


2 


.185 


it 


0. 


,007 


0, 


,762 


+0 
-0, 


098 
,068 


0, 


1 cn+O. 

.160+0 


,053 
,053 


0, 


,093 


+0 
-0 


,015 
,(315 


10, 


.48 


it 


0, 


.34 


NGC0720 


E 


2 


.381 


it 


0, 


Oil 


0, 


,653 


+0. 
-0, 


326 
,174 


0, 


AAO + 0- 


,114 
,083 


0, 


,249 


+ (3 
-0 


,056 
,025 


11 


.16 


it 





.27 


NGC0821 


E 


2 


.281 


± 


0, 


,005 


0, 


,950 


+0. 
-0 


068 
,038 


0, 


01 0+0! 
,210_J5 


,038 
,023 


0, 


137 


+0 
-0 


,015 
,005 


10, 


.85 


it 





.28 


NGC0936 


SO 


2, 


.258 


± 


0, 


,026 


1, 


,110 


+0. 
-0, 


098 
,098 


-0, 


.084Ip 


,053 
.053 


0, 


,121 


+0 
-0 


,035 
,015 


10, 


.27 


it 





.42 


NGC1316 


SO 


2, 


.344 


it 


0, 


,024 


0, 


,510 


+0. 
-0, 


.038 
,053 


0, 


,260tQ; 


,053 
.038 


0. 


,109 


+(3 
-0 


,025 
,015 


11 


.41 


± 


0, 


.22 


NGC1336 


E 


1, 


.982 


it 


0, 


,024 


1, 


,073 


+0. 
-0, 


,114 
,068 


-0, 


.3673q 


.053 
.053 


0. 


,122 


+0 
-0 


,035 
,015 


10, 


.02 


it 


0, 


.31 


NGC1339 


E 


2, 


.199 


it 


0, 


,024 


1, 


,146 


-t-O. 

-0, 


.098 
,129 


-0 


.079+0 


.053 
.053 


0. 


,173 


+(3 
-0 


,035 
.015 


10, 


.19 


it 





.40 


NGC1351 


E 


2 


.196 


± 


0, 


,024 


1, 


,129 


+0 
-0 


098 
.083 


— 


1 1 -7+0, 

.117+0 


,053 
,038 


0. 


,140 


+0, 
-0, 


.025 
.016 


10, 


.50 


± 


0, 


.31 


NGC1373 


E 


1 


.875 


it 


0, 


,024 


1, 


,010 


+0 
-0, 


(J83 
,053 


-0, 


01 0+0 
.313_(J 


,068 
,083 


0. 


,111 


+0 

-0, 


.035 
.016 


9, 


.45 


it 


0, 


.45 


NGC1374 


E 


2, 


.267 


it 


0. 


,024 


0, 


,931 


+0. 
-0, 


083 
,068 


0, 


1 1 /+0. 

■114lo 


068 
,114 


0. 


,149 


+0, 
-0 


.025 
.016 


10, 


.59 


it 


0, 


.31 


NGC1375 


SO 


1, 


.748 


it 


0. 


,072 


0, 


,323 


+0. 
-0, 


.023 
,023 


0, 


.068+0 


,053 
,068 


— 0. 


,045 


+0, 
-0 


.035 
.025 


9, 


.34 


it 


0, 


.40 


NGC1379 


E 


2 


.114 


it 


0, 


,024 


0, 


,913 


+0. 
-0, 


.083 
,053 


— 


.020+g 


,053 
.053 


0. 


,154 


+ 
-0 


,025 
,015 


10, 


.26 


it 





.32 


NGC1380 


SO 


2 


.340 


± 


0, 


,024 


1, 


,008 


+0. 
-0, 


114 

,083 


0, 


1 C-T + O. 

.157^0 


,098 
.068 


0. 


,127 


+0 

-0 


,035 
,015 


10, 


.85 


± 


0, 


.28 


NGC1380A 


SO 


1 


.740 


it 


0, 


,072 


0, 


,476 


+0. 
-0, 


.235 
,083 


-0, 


.012+° 


,114 
,098 


-0. 


,041 


+0 
-0 


,056 
.035 


9, 


.21 


it 





.40 


NGC1381 


so 


2, 


.185 


it 


0, 


,024 


0, 


,923 


+0. 

-0, 


.068 
,038 


0, 


1 A +0 

.104+JJ 


(353 
,038 


0. 


,092 


+0 
-0 


.025 
.005 


10 


.12 


it 





.39 


NGC1399 


E 


2, 


.574 


± 


0, 


,024 


1 


,026 


+0 
-0 


()83 
.053 


0, 


.346+0 


(368 
,038 


0. 


212 


+0 

-0 


,025 
,015 


11 


.91 


± 





.21 


NGC1404 


E 


2 


.415 


± 


0. 


,024 





,994 


+0 
-0 


()fi8 
.053 


0, 


.187+0 


068 
,038 


0. 


128 


+0, 
-0 


,025 
,005 


10, 


.89 


± 





.32 


NGC1419 


E 


2 


.068 


± 


0, 


,024 


1, 


,308 


+0 
-0 


220 
.220 


-0, 


.506+0 


,098 
,189 


0. 


173 


+0 

-0 


,045 
,025 


9, 


.74 


± 


0, 


.45 


NGC1427 


E 


2 


.243 


it 


0, 


,024 


0, 


,981 


+0. 
-0, 


.038 
,023 


— 0, 


.062t° 


,023 
.023 


0. 


,104 


+0, 
-0, 


,015 
.006 


10, 


.72 


it 


0, 


.29 


NGC1453 


E 


2, 


.475 


it 


0, 


,005 


1, 


,008 


-t-O, 

-0, 


,053 
,053 


0, 


.-1-0 

■294Iq 


.038 
.023 


0. 


,182 


+0 
-0 


.015 
.006 


11 


.45 


it 





.39 


NGC1461 


SO 


2 


.294 


± 


0, 


,011 


0, 


,864 


+0. 
-0, 


.098 
,098 


0, 


-1 r r + O! 

.155+U 


,083 
,053 


0. 


134 


+0, 
-0 


.025 
.015 


9, 


.88 


± 


0, 


.67 


NGC1600 


E 


2 


.508 


± 


0, 


,009 


0, 


,822 


+0. 
-0 


.189 
,114 


0, 


.409lo 


(353 
,023 


0. 


224 


+0, 
-0 


,025 
,015 


11 


.84 


± 


0, 


.36 


NGC1700 


E 


2 


.356 


± 


0, 


,005 


0, 


,519 


+0 
-0 


()38 
.008 


0, 


nor + O 

.335+0, 


(323 
,023 


0. 


114 


+ (3 

-0 


,015 
,005 


11 


.09 


it 


0, 


.32 


NGC2300 


E 


2 


.418 


it 


0. 


,006 


0, 


,719 


+0 
-0, 


053 
,144 


0, 


.404+0 


,053 
,023 


0. 


,212 


+0, 
-0 


,015 
,005 


11 


.13 


it 





.38 


NGC2560 


SO 


2, 


.303 


it 


0, 


,006 


0, 


,810 


+0. 
-0, 


159 

,174 


0, 


.226to 


,068 
,053 


0. 


,096 


+0 
-0 


,025 
,015 


10, 


.26 


it 





.94 


NGC2778 


E 


2 


.198 


± 


0, 


,007 


0, 


,897 


+0. 
-0, 


.098 
,174 


0, 


.243+0 


,053 
,053 


0. 


,166 


+0, 
-0 


,025 
,005 


10, 


.39 


± 


0, 


.35 


NGC3115 


SO 


2 


.378 


it 


0, 


,053 


0, 


,574 


+0. 
-0, 


.159 
,068 


0, 


.384^0 


053 
,053 


0. 


,095 


+0, 
-0 


,025 
,025 


10, 


.56 


it 


0, 


.32 


NGC3377 


E 


1, 


.981 


it 


0. 


,009 


0, 


,653 


+0. 
-0, 


.068 
,098 


0, 


.038+0 


(338 
,038 


0. 


,131 


+ (3 
-0 


,025 
,005 


9, 


.90 


it 


0, 


.28 


NGC3379 


E 


2, 


.291 


it 


0. 


,005 


0, 


,969 


+0. 
-0 


.038 
,023 


0, 


.i59t°: 


,038 
,023 


0. 


,168 


+0 
-0 


,015 
,(3(35 


10, 


.50 


it 


0, 


.28 


NGC3384 


SO 


2, 


.140 


it 


0, 


,031 


0, 


,522 


+0. 
-0, 


083 
.053 


0, 


.332+0 


,053 
,038 


0. 


,045 


+0 
-0 


,025 
,015 


9, 


.70 


it 





.44 


NGC3412 


SO 


2, 


.013 


± 


0, 


,034 


0, 


,632 


+0 
-0, 


098 
,159 


0, 


.029+0 


,083 
,083 


-0. 


Oil 


+0 

-0 


,045 
,035 


9, 


.34 


± 


0, 


.49 


NGC3585 


so 


2, 


.307 


it 


0, 


,122 


0, 


,710 


+0, 
— 


129 
.159 


0, 


.247+?; 

— 0. 


,068 
,053 


0. 


,060 


+0, 
— 


,035 
.015 


11, 


.03 


it 


0, 


.35 


NGC3607 


so 


2, 


.332 


it 


0. 


,101 


0, 


,720 


+0, 
-0, 


189 
,174 


0, 


.285+°' 


(368 
,068 


0. 


,098 


+0, 
-0 


.035 
.026 


11, 


.07 


it 


0, 


.33 


NGC3608 


E 


2, 


.259 


it 


0. 


,006 


0, 


,930 


+0, 
-0, 


,083 
,083 


0, 


.224+°' 


,053 
,038 


0. 


,136 


+0, 
-0 


.015 
.005 


10 


.78 


it 


0, 


.28 


NGC3818 


E 


2, 


.246 


it 


0, 


,007 


0, 


,763 


+0, 
-0, 


,144 
,068 


0, 


.361+0' 


,038 
,023 


0. 


,186 


+ 0, 
-0 


,015 
,005 


10, 


.73 


it 





.35 


NGC3941 


SO 


2, 


.117 


± 


0, 


,036 


0, 


,422 


+0, 
-0, 


,083 
,098 


0, 


.31lt°: 


,083 
,098 


0. 


,087 


+0, 
-0 


,045 
,035 


9, 


.71 


± 


0, 


.42 


NGC4026 


so 


2, 


.258 


it 


0, 


,029 


0, 


,542 


+0, 
-0, 


159 
,083 


0, 


.1531°: 


,098 
,083 


0. 


,088 


+0, 
-0 


,056 
.035 


10, 


.08 


it 





.41 


NGC4036 


so 


2 


.220 


it 


0, 


,039 


0, 


,574 


+0, 
-0 


114 

,129 


0, 


.4181°: 


083 
,038 


0. 


,253 


+0, 
-0 


.025 
.015 


10 


.33 


it 





.45 


NGC4111 


so 


2, 


.127 


± 


0, 


,019 


0, 


,464 


+0, 
-0 


023 
.008 


0, 


.074+°' 


,068 
,083 


0. 


109 


+0, 
-0 


,045 
,035 


9, 


.56 


it 


0, 


.56 


NGC4251 


so 


2, 


.072 


it 


0, 


,074 


0, 


,514 


+ 0, 
-0 


()3S 
.038 


0, 


■1731°: 


,053 
,038 


0. 


,040 


+ 0, 
-0 


,015 
,(315 


10, 


.09 


± 


0, 


.38 


NGC4261 


E 


2, 


.479 


± 


0, 


,005 


1, 


,136 


+ 0, 
-0 


053 
.053 


0, 


.2841°: 


,038 
,023 


0. 


175 


+ (3 

-0 


,015 
,(3(35 


11 


.40 


± 





.28 


NGC4697 


E 


2, 


.211 


it 


0, 


,009 


0, 


,924 


+0, 
-0, 


,068 
,053 


0, 


.027+°' 


,053 
,038 


0. 


,090 


+0, 
-0, 


,025 
.006 


10, 


.72 


it 


0, 


.23 


NGC5638 


E 


2 


.186 


it 


0, 


,006 


0, 


,946 


+0, 

-0, 


,068 
,038 


0, 


.1851°: 


,063 
,023 


0. 


,162 


+0, 
-0 


.015 
.005 


10 


.68 


it 





.29 


NGC5812 


E 


2, 


.300 


± 


0, 


,005 


0, 


,602 


+0, 

-0, 


,235 
,114 


0, 


.403+°' 


,038 
,098 


0. 


188 


+0, 
-0 


.005 
.025 


10, 


.72 


± 


0, 


.33 


NGC5813 


E 


2, 


.327 


± 


0, 


,007 


1, 


,127 


+0, 
-0, 


,098 
,083 


0, 


.1001°; 


,053 
,038 


0. 


,197 


+0, 
-0 


.015 
.006 


11, 


.21 


it 


0, 


.26 
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Table Bl. Continued. 



Galaxy Type Log((r) log(Age) [Z/H] [E/Fe] log(Mdy„) 



NGC5831 


E 


2, 


.206 


± 


0, 


,004 


0, 


,560 


+0, 
-0 


,023 
.023 


0, 


,360l"; 


023 
OOS 


0, 


,i35i;;: 


015 
,005 


10, 


,63 


± 


0.31 


NGC5846 


E 


2, 


.354 




0, 


,005 


0, 


,980 


+0, 
-0 


,068 
.038 


0, 


,296+0 


038 
023 


0, 


,2i5i'>: 


015 
(105 


11, 


,38 


± 


0.23 


NGC5866 


SO 


2, 


,143 


± 


0, 


116 


0, 


,413 


+0, 
-0, 


,098 
.098 


0, 


,2151°: 


114 
,114 


0, 


,089+0 


056 
045 


10, 


,39 


± 


0.36 


NGC6127 


E 


2, 


.393 


it 


0, 


,007 


0, 


,981 


+0, 
-0 


,053 
.038 


0, 


,266+0 


,038 
,023 


0, 


,217+°' 


015 
,005 


11, 


,25 


± 


0.41 


NGC6702 


E 


2, 


.243 


± 


0. 


,005 


0. 


,325 


+0, 
-0 


,038 
.008 


0, 


,439+0' 


,038 
,023 


0, 


,102+0' 


,015 
,015 


11. 


,00 


± 


0.30 


NGC6703 


E 


2. 


.258 


± 


0. 


,006 


0. 


,790 


+0, 
-0 


,083 
,068 


0, 


,180+°' 


053 
()23 


0, 


,110+°' 


015 
005 


10. 


,68 


± 


0.32 


NGC7052 


E 


2, 


.464 




0, 


,006 


1, 


,218 


+0, 
-0 


,068 
.068 


0, 




038 
,038 


0, 


,213+°' 


015 
,005 


11, 


,54 


± 


0.38 


NGC7454 


E 


2, 


,028 


± 


0, 


Oil 


0, 


,692 


+0, 
-0 


,053 
.144 


-0, 


■1281°: 


,053 
,038 


0, 


,031+°' 


,025 
,015 


10, 


,20 


± 


0.31 


NGC7562 


E 


2, 


,391 


± 


0, 


,003 


0, 


,924 


+0, 
-0 


,023 
,023 


0, 


,2031°; 


,038 
008 


0, 


,144+°' 


015 
()()5 


11, 


,29 




0.35 


NGC7619 


E 


2, 


,498 


± 


0, 


,004 


1, 


,051 


+ 0, 
-0 


,053 
,023 


0, 


,320+0 


023 
()23 


0, 


,i73i;;: 


005 
()()5 


11, 


,58 


± 


0.30 


NGC7626 


E 


2, 


,414 


± 


0, 


,004 


0, 


,963 


+ 0, 
-0 


,038 
,023 


0, 


■3361°: 


023 
008 


0, 


,2i8+;;' 


005 
005 


11, 


,36 


± 


0.,37 


NGC7785 


E 


2, 


,378 


± 


0, 


,005 


0, 


,911 


+0, 
-0 


,068 
.038 


0, 


,2661°: 


038 
,008 


0, 


,1511°: 


,015 
,005 


11, 


,19 


± 


0.39 



